A Hartree-Fock Study of Charge Redistribution in a 2D Mesoscopic Structure 
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In tliis paper, we investigate the ground state of two-dimensional disordered cylinders which contain 
spinless, interacting electrons using the Hartree-Fock approximation. Calculations of the deviation 
of the polarization from uniformity reveal a tendency of the charge to rearrange towards the ends of 
the system. The presence of disorder results in fluctuations of the deviation around its mean value, 
which are more pronounced when the disorder strength is of the order of the interaction between 
the electrons. 



The existence of persistent currents in normal meso- 
scopic rings threaded by a magnetic flux j|] has stimu- 
lated a great deal of experimental and theoretical work, 
lH- More recently, there have been studies of the 

ground state of one-dimensional rings containing spin- 
less fermions ( [Q- [^) which take into account both 
electron-electron interactions and disorder. Strongly dis- 
ordered rings with short-range interactions and a half- 
filled band exhibit a reorganization of the charge in the 
ground state. The charge density changes from an in- 
homogeneous configuration due to the presence of the 
strong disorder to a periodic array of charges as a result 
of the interactions. 

In the present paper, we study the ground state of 
two-dimensional cylinders which contain spinless, inter- 
acting electrons via Coulomb interactions, by solving self- 
consistently the Hartree-Fock equation. The validity of 
the method has already been established ||l^ by studying 
the ground state of one-dimensional disordered rings and 
finding good agreement with the exact calculations [p^ . 
Our aim is to examine reorganization of the charge in 
the ground state, and in particular, the manner in which 
charge tends to polarize towards the ends of the cylinder, 
thereby producing spontaneous dipole or quadrapole mo- 
ments. 

The system under investigation is a two-dimensional 
cylinder, formed by a tight-binding lattice with M sites 
in the longitudinal direction (x) and L sites in the trans- 
verse direction (y). We take periodic boundary condi- 
tions along the longitudinal direction and free boundary 
conditions along the transverse. The cylinder contains N 
spinless electrons described by the Hamiltonian 
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where each site I has coordinates I — {x,y). The opera- 
tors cj, ci create and destroy a particle at site /, respec- 
tively. Vik is the hopping element between different sites. 
In the following, we will restrict ourselves only to nearest- 
neighbour hopping elements of strength Vik = —V . ei is 
the on-site energy which is equal to e; = 4:V -\-ri, where r; 
are random numbers uniformly distributed over the range 
[-W/2,+W/2]. W is the strength of the disorder, and 
in the clean case W = 0. Uik is the interaction between 
the particles, which has been taken to be long range. 
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where r„ is the position of the nth particle. The Hartree- 
Fock equation which corresponds to the Hamiltonian (|l|) 
is of the form 
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where is the amplitude of the nth single-particle 

wavefunction on site I, and En is the corresponding 
single-particle energy. The third and fourth terms are 
the direct and exchange potentials, respectively. 

The probability of finding an electron on site (l) is 
'^x,y = J2n=i ^^^^ heuce the mean number of 

electrons per unit length in the transverse direction is 



M 



piy) = E ' 



(4) 



1 



We then define 
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where p{y)Ay is the probability of finding an electron in 
the interval [y,y + Ay]. A measure of the deviation of 
the charge from uniformity is the standard deviation of 

p{y) 
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where < y" >= ^ y"p{y)Ay. In the case of a clean 
system {W = Q) the electron density for a uniform dis- 
tribution must be p{y) = N/L, where L is the system 
size in the transverse direction (number of chains). This 
gives a^/L^ ~ j^. On the other hand, for the extreme 
case in which half the charge is at L/2 and the other 
half at -L/2, p{y) = fS{y + f ) + f <5(y - f ) which 
gives / = Thus, the length-normalized deviation 
of p{y) can take values in the interval [1/12,1/4], the 
minimum and the maximum values corresponding to a 
uniform distribution and a maximally "polarized" one, 
respectively. 

In figure ^ we present the normalized deviation with 
respect to the number of particles N, for different values 
of U, and V — 1. The system has no disorder and its size 
is L — AI — 10. As we increase the interactions between 
the particles, for small N, cr^/L^ aproaches its maximum 
value, which indicates total "polarization". This is the 
expected result since we are approaching the electrostatic 
limit. When the system contains more electrons there 
is still an increase in cr^/L^ as U increases, but not all 
electrons move towards the ends of the cylinder. We 
have confirmed that this again approaches the expected 
minimum electrostatic energy configuration. 

In figures ^ and ^ we show for a system with non- 
zero disorder. In figure ^ we plot the ensemble average 
of versus U, for disorder strength W = 2 and iV = 10 
along with for three individual samples. Individual 
samples show small fluctuations around the mean value 
of (7^ indicating that the charge tends to separate. In- 
creasing the disorder, as in figure ^ where W — A, one 
can see that the fiuctuations increase with disorder, even 
though the average (dots) remains practically unaffected 
by it. However, the behaviour of the first moment < y > 
for individual samples is more pronounced than that for 
a^. In figures |^ and || we have plotted < y > for two 
different strengths of disorder, W = 2 and W = A, re- 
spectively. In the inserts of ^ and|^, the ensemble average 
of the absolute value of the first moment is shown. 

Figures ^, ||, § and || illustrate behaviour which we ex- 
pect to be typical of small systems with free-end bound- 
ary conditions in at least one direction. Namely that with 
increasing U the dipole moment < y > induced by ran- 
dom potential fluctuations (ie. W ^ 0) tends to zero, and 
that the charge distribution becomes sharply peaked at 



the ends of the samples. A key feature revealed by these 
flgures is that at intermediate values of U (of order W), 
both and < y > exhibit large sample-to-sample fluc- 
tuations about their means and large fluctuations with 
increasing U, associated with charge redistribution of the 
ground state. 

Finally, in figures ||, |^ and || we present cr^ , < y > and 
the current /, respectively, as a function of the phase 
(f> — 27r(f>/<I>o, where is the magnetic flux threading the 
cylinder, and <I>o is the flux quantum. The current is 
equal to 
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where Eg is the Hartree-Fock ground state energy. The 
results show the behaviour of one disordered cylinder, for 
the case of = 2 and for one value of the interaction 
strength ([/ ~ 4) such that it is of the order of the dis- 
order strength. A common feature in these three figures 
is that all quantities are symmetrical around the value 
of the phase cj) which corresponds to a magnetic field of 
half a flux quantum. Such charge fluctuations could pos- 
sibly be detected experimentally by placing a SET in the 
vicinity of a sample, which couples to the electric field 
generated by such a non-uniform charge distribution. 

In this paper, we have made a Hartree-Fock study of 
the ground state of two-dimensional cylinders which con- 
tain spinless, interacting electrons. We found that the 
charge in the ground state shows a separation towards 
the ends of the cylinder, which is reflected in the polar- 
ization of the system. The polarization shows fluctua- 
tions around its mean value in the presence of disorder. 
These fiuctuations are stronger in the first moment of the 
charge disrtibution when the disorder, bandwidth and 
interaction between the electrons are of the same order, 
refiecting the competition between Mott and Anderson 
localization ( p8[). 
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FIG. 1. Length-normalized deviation of p{y) with respect 
to the number of particles A'^ for a clean cylinder L = M = 10, 
for different values of the interaction strength U. 
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FIG. 2. Length-normalized deviation of p{y) with respect 

to the interaction for disordered cylinders of size L = M = 10 
and A*' = 10 electrons, with W — 2. The dots represent 
the ensemble average of obtained from 200 samples. The 
different lines represent individual disorder realizations. 
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FIG. 3. Length-normalized deviation of p{y) with respect 

to the interaction for disordered cylinders of size L — M = 10 
and A'^ = 10 electrons, with W — 4. The dots represent 
the ensemble average of obtained from 200 samples. The 
different lines represent individual disorder realizations. 
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FIG. 4. First moment of versus U for disordered cylin- 
ders of size L = M = 10 and iV = 10 electrons, with W = 2. 

The dots represent the ensemble average. The different lines 
represent individual disorder realizations. The insert shows 
the average of | < j/ > | versus U. 
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FIG. 5. First moment oi p(y) versus U for disordered cylin- 
ders of size L = M = IQ and A'' = 10 electrons, with W = A. 
The dots represent the ensemble average. The different lines 
represent individual disorder realizations. The insert shows 
the average of | < y > | versus U . 
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FIG. 8. Current versus magnetic field for an individual dis- 
order realization. The cylinder is of size L = M = 10 and 
contains N — 10 electrons, with W = 2. The results have 
been obtained for U = 4. 
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FIG. 6. Length-normalized deviation of for an individ- 
ual disorder realization, as a function of the magnetic field. 
The cylinder is of size L = M = 10 and contains iV = 10 
electrons, with W = 2. The results have been obtained for 
f7 = 4. 
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FIG. 7. First moment of p{y) for an individual disorder 
realization, with respect the magnetic field. The cylinder is 
of size L = M = 10 and contains A'^ = 10 electrons, with 
W = 2. The results have been obtained for U = 4. 
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